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Using Gaussian null coordinates, the near-horizon spacetime metric is pertubatively constructed 
for both isolated and dynamical trapping horizons (equivalently marginally outer trapped tubes). 
Spacetime is allowed to be of arbitrary dimension and the formalism accommodates both general rel- 
ativity and more general field equations. The formalism is demonstrated for two applications. First, 
spacetime is considered near an isolated horizon and the construction is both checked against the 
Kerr-Newman solution and compared to the well-known near-horizon limit for stationary extremal 
black hole spacetimes. Second, spacetime is examined in the vicinity of a slowly evolving horizon 
and it is demonstrated that there is always an event horizon candidate in this region. The geome- 
try and other properties of this null surface match those of the slowly evolving horizon to leading 
order and in this approximation the candidate evolves in a locally determined way. This general- 
izes known results for Vaidya as well as certain spacetimes known from studies of the fluid-gravity 
correspondence. 



I. INTRODUCTION 

There is a large literature on geometrically denned black holes and their horizons. Trapped surfaces and apparent 
horizons were defined close to fifty years ago £Q, but in the last two decades interest in both the mathematical and 
physical properties of these objects has increased. Almost twenty years have passed since Hayward's original definition 
of trapping horizons [2 while isolated and dynamical horizons have been studied for over a decade 3 , 4 . In that time 
these geometric horizons have been widely studied in mathematical relativity (see for example [5Tl37|) but have also 
found applications in numerical relativity (examples include |38H47j ) loop quantum gravity (for example [48H50] ) , and 
the fluid-gravity duality regime of the AdS-CFT correspondence [5D - f54] . 

Very recently work has begun to study spacetime close to the horizon[51 [SS]. This paper presents the details of 
results first announced in |55j . developing the necessary mathematics for stepping off of the horizon and studying 
the near-horizon spacetime. We expect this work to find many uses, but here we concentrate on the mathematical 
formalism followed by just two applications. In the first we construct the spacetime around an extremal isolated 
horizon and show that at leading order it takes the familiar near-horizon form found in works such as |56l 157) . In 
the second we construct spacetime around a general slowly evolving horizon and demonstrate the existence of a null 
surface that hugs the horizon. This candidate event horizon has previously been found for specific spacetimes including 
Vaidya [TTJ [T2J [M] and several black-brane spacetimes that show up in the fluid-gravity correspondence [5TM54) . 

Our construction applies to spacetimes of arbitrary dimensions with horizons of arbitrary signature. However, 
in essence it is a generalization of the standard (3 + l)-formulation of general relativity. In that formalism the 
Einstein equations are decomposed into constraints on the intrinsic and extrinsic geometry of a three-surface, along 
with evolution equations that determine how that geometry changes as the surface is propagated forwards in time. 
Similarly we will derive a set of constraint equations for the allowed geometry of a horizon along with evolution 
equations that describe how that geometry changes as one moves away from the horizon. These are used to construct 
spacetime in a neighbourhood of the horizon. 

Mathematically we work with (n+l)-dimensional spacetimes and consider n-dimensional hypersurfaces which can 
be foliated into (n — l)-dimensional spacclike surfaces. This mirrors the structure of dynamical trapping horizons 
which come with a unique foliation [32] as well as the marginally outer trapped tubes found in numerical relativity 
(which are built from (n— l)-dimensional apparent horizons found on successive time slices). How the geometry of 
the (n— l)-leaves changes under deformations was studied in detail in |22j and we import many results from that 
paper. However, that work must be supplemented in order to understand the full geometry of horizons and other 
nearby n-dimensional surfaces. We also need to understand how the slices fit together and how that fit changes under 
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deformations. 

After studying the geometry of deformations in some generality we specialize to find the near-horizon spacetime 
metric in Gaussian null coordinates constructed off of the horizon. Given a horizon H and its foliation S v , we take the 
(inward) null normals to the S v and consider the null geodesies that are tangent to those normals at H . Coordinates 
are then constructed taking the affine parameter p along the geodesies as the off-horizon coordinate and Lie-dragging 
the on-horizon coordinates to the level surfaces of p. The components of the metric in these coordinates can be 
calculated (perturbatively) by considering how the geometry changes if H is deformed by the d/dp vector field. With 
the near-horizon geometry in hand, we move on to the aforementioned applications. 

The paper is organized as follows. In Section [II] we briefly review the (n + l)-formulation of general relativity. 
The intent of this section is mainly to demonst rate the processes that will be used in later sections but in a setting 
that is more familiar to most readers. Section III reviews the geometry of (n— I )-dimensional surfaces and the n- 
dimensional hypersurfaces that can be built from them and then applies that work to perturbatively reconstruct 
spacetime near such hypersurfaces. Next, Section IV reviews definitions of the various types of geometric and causal 
horizons. Section [V] applies the definition of an isolated horizon to reconstruct the spacetime around an extremal 
isolated horizon while Section |VI| works from slowly evolving horizons to demonstrate the existence of a horizon- 
hugging event horizon candidate. As a cross-check it also compares this general result with analogous ones known 



for particular spacetimes. Section [VII| summarizes this work and considers some future applications. Finally, working 
with Kerr-Newman spacetimes, Appendix [A] demonstrates both the construction of Gaussian null coordinates around 
an isolated horizon as well as how the near-horizon metric may be reconstructed using horizon data (plus a bit more) . 

Note: Another programme very similar to this one has been independently initiated by Krishnan [5]. As in our 
approach he uses Gaussian normal coordinates and rebuilds spacetime metrics (to second order) from the deformations 
of quantities such as the expansions and shears. However the appearance of the mathematics is quite different as 
Krishnan's paper is written in the language of the Newman-Penrose formalism. It specialized to isolated horizons in 
four-dimensional spacetimes, though at the same time it explicitly incorporates electromagnetic fields. Which of the 
formalisms is more useful will depend on the details of a particular application. 



Notation 



Throughout we assume an (n+l)-dimensional spacetime {M,g a f,,'V a ) and study embedded surfaces of dimensions n 
and (n— 1). We use lower-case early-alphabet latin letters {a, b, c, . . . , g} as abstract indices on the full spacetime but 
switch to greek letters {a, /3,7. . . } when working with a coordinate chart. Similarly lower-case mid-alphabet latin 
letters {h,i,j, . . .p} are used as abstract indices for tensors in n-dimensional surfaces while sans-serif versions of the 
same letters {h, i,j, . . . p} are used for coordinates and tensor components relative to coordinates. Finally upper-case 
latin letters {A, B,C . . .} are used as abstract indices on (n— l)-dimensional surfaces while their sans-serif versions 
{A, B, C . . . } are used to indicate quantities written in terms of coordinates. 

The pull-back operator between surfaces will always be written as an e with indices indicating which spaces it 
operates between. Thus the induced metrics on n- and in— l)-dimensional surfaces are respectively 

qij = efe^gat and q AB = e a A e b B g a b , (1) 

while tangent vectors to those same surfaces would push-forward to TM via 

N a = e?V l and V a = e a A V A . (2) 

If we switch to coordinate charts so that (sections of) n- and (n — l)-dimcnsional surfaces are parameterized by 
functions 

x a = X a (y') and x a = Z a {9 k ) , (3) 

then 

e» - — and e" — U) 

e < dy< and " A ~ ■ (4) 

Fixing i or A, these are also the components of the coordinate tangent vectors d/dy' and d/d8 A (after they have been 
pushed-forward to TM). 

We follow the sign conventions of Wald [SI] for such things are Riemann and extrinsic curvatures. 
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II. GEOMETRY OF SPACELIKE HYPERSURFACES 



To establish basic ideas about hypersurfaces and how their geometry changes under deformations we begin with a 
quick review of the (n+1) -formulation of general relativity. For more details see standard texts such as [59H6T] . 



A. Basic geometry of an n-dimensional spacelike hypersurface 

Let (E, qij, be a spacelike n-dimcnsional surface embedded in an (rt+l)-dimensional spacetime (M, g ab , V Q ). As 
noted in the preamble, the induced metric on E is 

lij = Cttfgab , (5) 

while the corresponding extrinsic curvature is 

Kij = e-e^Wah , (6) 

where f a is the future-oriented unit normal to E. 

As for the elementary differential geometry of surfaces in Euclidean IR 3 , the intrinsic metric and extrinsic curvature 
are not independent but instead are related to each other as well as the curvature of the ambient spacetime through 
the Gauss-Codazzi equations. From the Gauss relation one can show that 

G ab r a f b = i (Re + K 2 - Kyi?*) , (7) 

while from the Codazzi relation 

e?G ab r b = DjK/ - D t K . (8) 

In these equations, i?g is the Ricci scalar for E, K — q l -> is the trace of its extrinsic curvature and G a b — r R>ab~\R-9ab 
is the Einstein tensor for M. Applying the Einstein equations 

G a b + A.g a6 = 8nT ab (9) 

to these, the Einstein tensor is replaced by terms involving the stress-energy tensor and cosmological constant. Then 
equations (J7| and ^ respectively become the Hamiltonian and momentum constraint equations. These conditions 
necessarily hold if E is surface embedded in a full solution of the Einstein equations. In particular, if (E, q^, K^) 
is viewed as in instantaneous configuration that will be time-evolved into a full spacetime (the viewpoint taken in 
numerical relativity) then the fields must satisfy these constraints to be valid initial data. 



B. Deforming an n-dimensional spacelike hypersurface 

Next, we consider how the induced metric and extrinsic curvature change if a hypersurface E D is deformed. This 
is easiest to understand if we switch to working with coordinate charts. Let {x a } be a set of (n+1) coordinates over 
(at least some region of) M and {y'} be a set of n coordinates so that an initial (portion of) E Q is parameterized by 
equations 

x a = X a (y [ ) . (10) 

Now a deformation may be defined by specifying a scalar field N and vector field N l over E Q . These are respectively 
called the lapse and shift and used to construct an rj-dimcnsional vector field 

T a = Nr a +N a = NT a + e?N\ (11) 

In turn this is used to (infinitesimally) generate a new surface defined by 

X a (y')^X a (y') + (At)T a (y'). (12) 

The mapping also identifies points on the hypersurfaces (essentially by Lie-dragging coordinates between surfaces) as 
shown in FIG. [T] Derivatives relative to T a are defined in terms of the usual limits. For example the deformation of 
the induced metric is 

STq -'i = i£b to ' (13) 



FIG. 1: Time evolution of a spacelike slice. The time-evolution vector T a deforms an initial surface E D into £At- It can be 
decomposed its parts perpendicular and parallel to £, hence defining a lapse function N and shift- vector field N a = efiV 1 . 



where the metrics on the two surfaces may be directly compared since they are expressed relative to the same (Lie- 
dragged) coordinate system. This is a covariant construction. 

Of course one doesn't use the limit definition to actually compute deformations. Instead calculations are based on 
the understanding that the deformations effectively extend the coordinate system {y'} off of S D and supplement it 
with another coordinate t such that 

r=|and„ = ^. (14) 

These are coordinate vectors and so their Lie brackets vanish. In particular 

£ T e?=0 (15) 

and in turn if follows that 

e?£ T T a = =► e«T b V b T a = -D t lnN + K tJ W . (16) 



Deformations are calculated by taking Lie derivatives with respect to T a with appropriate applications of ( 16 ) to 
convert derivatives normal to £ into derivatives tangent to it. The derivative of the induced metric can be calculated 
without this condition 

Srqa = f$e)C T g ab = 2NKy + C N q l3 (17) 
however the deformation of is a little more complicated. Starting from 

8 T K l3 =e c le)CT{^ah) (18) 
a certain amount of algebra along with an application of the Gauss relation gives: 



S T K l3 = DiDjN - N(R ZJ - 2K ik K j k + KKy) + C N K VJ + e$e) - -q tJ g ab ) G ab , (19) 



where Rij is the rt-dimensional Ricci tensor for the surface. 



Of course, the best known example of a deformation is the time evolution of initial data in the (3 + l)-formulation 
of general relativity. Given an initial surface S, the Einstein equations are equivalent to the constraints Q and ^ 
along with the time evolution of the extrinsic curvature (19). Though the details will differ, this is the perspective 
that we will adopt when studying the spacetime around a horizon - we will take the specification of the horizon as 
initial data and then try to evolve that data to construct spacetime in a neighbourhood of that horizon. 



C. Spacetime near a spacelike hypersurface 

We can also use this knowledge of deformations to perturbatively construct spacetime close to a spacelike surface. 
Working in geodesic normal coordinates based on E„ the spacetime metric takes the form 

ds 2 = -dr 2 + h li dy ] dy i , (20) 
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where r measures proper time along the geodesies and the hy s are the components of the spacelike n-dimensional 
metric on surfaces E r of constant r. It is straight-forward to expand this as a Taylor series (in r) around the initial 
data on E D . To second order we have 

ds 2 a -dr 2 + (tfj + Tf^jj + y ^ dy'dyi , (21) 

where g;j is the induced metric on E Q , 

^ = 5 f ^ = 2K Vs (22) 

and 

= Sf (5f*j) = -2(i?ij - 2if ik ^ k + iOC u ) + 2 (e?ef - \%9 af> \ G a p , (23) 

since geodesic coordinates have N — 1 and vanishing shift vector. 

Initial data on a spacelike surface fully determines the spacetime in its future domain of dependence |59l 16 1] . Our 
Taylor expansion is consistent with this result. Momentarily restricting attention to vacuum spacetimes that are 
solutions of the Einstein equations (so that G a p — 0) , knowledge of qij and 2£y on E Q allows us to calculate all orders 
of derivatives of q\j . The first derivative of induced metric is determined by the extrinsic curvature while the derivative 
of the extrinsic curvature is determined by the induced metric and extrinsic curvature. This closes the circle and so, 
based on the initial data, we can construct the Taylor expansion to all orders. 

III. GEOMETRY OF HYPERSURFACES OF ARBITRARY SIGNATURE 

The previous section has reviewed several important ideas. First, a surface is characterized by its intrinsic and 
extrinsic geometry and in general these are not independent: they are linked to each other and the curvature tensor of 
the full spacetime through the Gauss-Codazzi equations. Second, we can calculate how the surface geometry changes 
if the surface is deformed by a vector field: computationally this amounts to taking Lie derivatives of surface quantities 
with the extra condition that the Lie derivative of pull-back operators vanishes. Third, given a coordinate system on 
the surface and a preferred vector field (in this case the tangent vector field to the congruence of timelike geodesies 
normal to S) we can construct a coordinate system in a neighbourhood of £ and use the deformation results to 
construct the spacetime metric relative to that coordinate system. The intrinsic metric and extrinsic curvature tensor 
of £ are good initial data (provided that they satisfy the constraints). Via the deformation/time evolution equations 
they fully determine the vacuum spacetime metric close to the horizon. 

These ideas form the foundation of the (3 + 1) formulation of general relativity where an initial instantaneous 
configuration is evolved into a full spacetime via evolution equations. Ultimately we wish to mirror this construction 
where the initial data is a geometric horizon H: given its intrinsic and extrinsic geometry we would like to construct 
spacetime in a neighbourhood of H . Unfortunately the standard formalism is not sufficient for our purposes. General 
horizons can be null and even when they are spacelike we will often be interested in the regime where they are "almost" 
null. Thus we do not wish to base our formalism around a timelike unit normal vector f a : it is not well-defined for 
truly null surfaces and is inconvenient for "almost null" surfaces. 

The standard formalism is also not sufficient in that we wish to work with spacetimes of arbitrary dimensions as 
well as allow for generalizations of the Einstein equations. However these are relatively minor issues. Switching the 
formalism to higher dimensions is straightforward while the Einstein equations are not actually a fundamental part 
of the formalism. Most of the formalism is general to any spacetime and the Einstein equations are only used to 
constrain the potential spacetimes under consideration. 

This section reformulates the standard formalism to accommodate our goals. As such we consider n-dimensional 
surfaces that can be defined as the smooth union of a set of spacelike (n—1) -dimensional surfaces S v : H = {U„5„} for 
some range of the surface label v. In order to ensure maximum generality and applicability, at this stage we do not 
impose any restrictions on the overall signature of H and neither do we assume that it has any particular geometric 
properties (for example we do not assume that it satisfies any constraint equations). 

A. Geometry of (n— l)-dimensional spacelike surfaces 

We begin by reviewing the geometry of the (n—1) -dimensional building blocks of H. As in the previous section we 
split this into a consideration of basic geometry followed by a study of deformations. This mathematics is well-known 
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and has been derived and re-derived many times. That said, our immediate reference (in which many more details 
can be found) is [22) . 

1. Basic geometry of (n— 1)- dimensional surfaces 

Let (S, qAB i d A ) be an (n — l)-dimensional spacelike surface embedded in a time-orientable (n+ l)-dimensional 
(M, g ab , Va)- The induced metric on S is the pull-back of the full spacetime metric: 

QAB = e A e b B g ab (24) 

and this determines the full intrinsic geometry of S including the metric compatible covariant derivative cIa and the 
Riemann tensor Rabcd- In the particularly important case where S is two-dimensional (n = 3) we have 

Rabcd = {qacQbd - QadQbc) , (25) 

for the two-dimensional Ricci scalar R. 

The normal space to S is two-dimensional. It is spanned by a pair of future-oriented null vectors £ a and n a and we 
assume that properties of 5* and the spacetime are such that these can respectively identified as outward and inward 
pointing. The direction of these vectors is fixed, but their scaling isn't. We remove one degree of freedom by requiring 
that they be cross-normalized so that 

£-n = -\ (26) 

which leaves a single degree of rescaling freedom: 

£ ^ eJ £ and n^e~ f n, (27) 
for an arbitrary function /. However, independent of that particular choice of scaling we have 

~ab a b ~AB ab , nab , ,a /jfr /oo\ 

q — e A e B q = g + £ n + n I . (28) 

If n = 3 then the induced area form on S relative to that of the full spacetime is 

e AB = e A e b B (e abfg £f n 3 ) . (29) 

Similarly for n = 4 and 5 the volume forms are 

cabc = e a A e B e c c {e abcfg £ f n 9 ) and Iabcd = e a A e B e c c e d D (e abcdfg £ f n 9 ) . (30) 

The generalization to even higher dimensions is obvious. 

Just as extrinsic geometry of E was determined by how t a varied over the surface, the extrinsic geometry of S 
can be understood by considering how these null normals vary along the surface. We have the extrinsic curvature 
analogues: 

k [ AB = e a A e b B VJ b and k^ ] B = q A q\V a n b (31) 
as well as the connection on the normal bundle: 

Co A = -e a A n b V a l h . (32) 

These measures of the intrinsic and extrinsic geometry are related to each other and the geometry of the ambient 
spacetime by analogues of the Gauss, Codazzi, and Ricci equations (see [55] for details). 

Note however that compared to the codimension-one spacelike case, there is an extra complication in dealing with 
these extrinsic curvature quantities. For spacelike codimension-one hypersurfaces, the extrinsic curvature is a uniquely 



defined geometric quantity. Here things are a little more complicated. Under the rescalings defined by (27 1 

k AB ~^ efk< AB i k AB ~^ e ~ fk AB an ^ &A ~t &A + d A In / . (33) 

and so they have a gauge dependence. For now we accept this uncertainty but when constructing the spacetime 



around H in Section III B we will fix the gauge. 

The traces of the extrinsic curvatures and their trace-free parts are important enough to have their own names. 
For a general element X a — a£ a — r )n a of the normal space we write 



k [ AB ] = e A e B V a X b = t^Qab + vffl , ( ^ ) 



where the trace 9(x) = Q AB k A B is the expansion and the trace-free <J A X B is the shear. 
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2. Deforming an (n— 1) surface 



Just as we studied how the geometry of £ is changed by a deforming vector field we can also examine the deformations 
of S. The process is essentially the same. Parameterize S in terms of a coordinate chart 8 A : x a = X a {9 A ). Then, for 
a deformation vector field X a : 



X a {6 A ) -> X a {6 A ) + eX a (6 A ) 



(35) 



defines a new surface S e by deforming 5 a coordinate distance e in the direction X a . It identifies points on the 
(n— l)-surfaces (essentially by Lie-dragging coordinates between surfaces) and so we have 



C x e\ = 0. 



(36) 



We can then examine how the geometry changes under these deformations. For our purposes, it will be sufficient to 
restrict attention to normal deformations so that 



X a = at - 7 n a , 

for some functions a and 7. Then it is straightforward to see that 

SxQab = e a A e b B C x g a b 
and (dropping the indices) 

= 0(x)£ = (<x9(t) - 7#(n))e- 



- -2k$=2(ak% 



(37) 



(38) 



(39) 



for the volume- form |74j . It is then obvious why we call the traces expansions. They tell us how the volume elements 
change while the shears are the part of the evolution that deforms S (but does not change its volume). 

As was the case for n-dimensional hypersurfaces, deformations of the extrinsic curvature quantities are more com- 
plicated. In addition to applications of the Gauss, Codazzi and Ricci relations we again need to convert off-horizon 
derivatives into ones tangent to the S. As in the previous section the key to this is ( 36 1 and this time the analogues 
to ( 16 1 are: 



X a V a £ b = -d a "/ + w a 7 + Kxlb and 
X a S7 a n b = -d a a + u a a - K X n b , 



(40) 



where d a f = e^d^f, £> a = e^a)^, k x = —X a n b \7 a £ b and the index-reversed = 3 a fc e s9 AS - Note that k x is a gauge 
dependent quantity whose value depends on how the scaling of the null vectors changes off the original surface. Under 
rescalings (27) of the null vectors: 



kx k x +C x (lnf) 



(41) 



In Section flH B| we will also gauge-fix this quantity but for now leave it undetermined. Because we have only considered 
normal deformations, there are no extrinsic curvature terms in (40) as compared to (16). 

For details of the deformation calculations see [221, here we will just list results. First deforming the extrinsic 

(42) 



curvature k^ B we find that 
Sxk AB 



.(<) 



-d A d B l + 2LU( A d B} -f + K X k AB 
+a (k%kf C - e a A e b e B e d C abcd 



1 



(n-1) 



+7 



\ Rab + 2 



1 r 



(i)k { AB +0 {n) 



AB 



1 



e a A e B TZ ab + d(A&B) ~ u A Cu B 



where C abc d is the (n+l)-dimensional Weyl tensor and Rab is the (n— l)-dimensional Ricci tensor. In the usual way 
parentheses indicate a symmetrization of indices so, for example, 



W(A^S)7 = 2 (^A^b7 + ^BdAl) 



(43) 



(n) 

The deformation of the extrinsic curvature k ab is 
£x*4"b = d A d B a + 2Q( A d B )a- Kxk^l 

ry ( h {n)C - p a n b p c n d r 

7 I K AC K B e A n e B n °a 



(44) 



bed 



(n-l) 



qABUcdrfn 



-a 



1 - 1 
~ 2 Rab + - 2 



2k< C(A k B) C + \ e A e B'^-ab + d {A L) B) + UJ A LU B j 



Given (42) one can easily find (44) by interchanging I and n in and then sending a — » — 7 and 7 — > — a. 

As noted, null expansions are particularly important in discussions of geometric horizons and so we also list the 
deformations of the traces of the extrinsic curvatures: 



>XV(t) 



Kx6(£) — d 7 + 2oj d A j 
-7 



(45) 



Ml 2 



d A u A - I + \n ab q ab 



a lib 



n ab £ a £ 



(n-l) W 



and 



-Kx6{n) + d^a + 2Q A d A a 



(46) 



Ml 2 



-A 



i? 1 



-ah 



'7 



V)f + K ab n a n b + 



1 



where d 2 = d A d/ 



AB (") 
a (n) AB 



^4, = u> Wi, i? is the (n— l)-dimensional Ricci scalar for S v , \\cr^\\ = o~f^ o y A ' B and ||cr(„) 

Finally the variation of the connection one-form is 

Sxua = -9(x)U A + d A K X - d B k A X g ] + ad A 6^) + jd A 6( n ) + e A TZ ab Xl . 



(47) 



Here Xf = al a + "/n a is normal to X a while 9x and k^jp are defined in the obvious way. If consulting [52] for details 
of this particular calculation, note that the last line equation (2.26) of that reference misses an overall factor of 1/2. 



B. Building an n-dimensional hypersurface from (n— l)-surfaces 

Next, we put these (n — l)-dimensional pieces together into an ?i-dimcnsional surface. As is now our standard 
procedure, we start by considering the intrinsic and extrinsic geometry and then move on to consider deformations of 
that geometry. 



1. Intrinsic geometry 

Let H be an n-dimensional hypersurface which is the union of a set of spacelike (n— l)-dimensional surfaces S v 
with v labelling the surfaces. Regardless of the signature of H, we can define a unique vector field V 1 on H that is 
1) normal to the S v , 2) tangent to H and 3) satisfies Cyv = 1. Then a general evolution vector field Y l that evolves 
leaves of the foliation into each other may be written as 

T l = V 1 + V , (48) 

where V 1 is tangent to the S v . The canonical example of an evolution vector field arises if we impose a foliation- 
compatible coordinate system z' = {v,9 A } on H. Then 

y = l (49) 

is an evolution vector field. For such a choice we have C y8 A = and so it follows that 



V A = -£ V 9 A . 



(50) 
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FIG. 2: This figure is similar in appearance to FIG. [T] and represents a similar situation. The bottom sheet H is foliated by 
(n — l)-dimensional surfaces S v (drawn as solid lines). On H, "V a is the evolution vector field evolving surfaces into each other 
from the left-hand to the right-hand sides and identifying points on those surfaces. The null normal n a is used to deform H 
into Ha p - The evolution vector field r Y a on Ha p will usually no longer be normal to the S v after the translation. 



In analogy with the evolution vector field in the (n+l)-formulation we will refer to V A as a shift vector field. 

Regardless of the value of the shift vector, the foliation parameter can be used to fix the scaling of the null vectors. 
We scale the null vectors so that 

v a = t _ Cn a > ( 51 ) 

for some expansion parameter C|75j. Note that if C > 0, H is spacelike, while if C < it is timelike and if C — it 
is null. Given this construction, the scaling freedom of the null vectors is restricted to the freedom to reparameterize 
the foliation labelling. For an alternative labelling v = v(v) we would have 

[dv] a = -jrHa =*■ V a = a(v)V a (52) 
a(v) 

where a(v) = is constant over each S v . Then 

T = at , h a = -n a and C = a 2 C . (53) 
a 

Under this restricted class of rescalings (with a constant over each individual S v ), uj a is invariant. 

Next consider the intrinsic and extrinsic geometry of the full H. Relative to {z' — (v, 6* A )} the intrinsic metric is 

dS 2 = q ij dz i dz j = 2Cdv 2 + 2q AB V A dvd6 B + q AB d9 A d6 B , (54) 

where qab is, as usual, the induced metric on the (n— l)-surfaces. Note that this coordinate form explicitly demonstrates 
how the sign of C determines the signature of H . As would be expected, the inverse metric g y is not well-defined if 
C — (that is when H is null). 

Shifting our attention to the extrinsic geometry, a future-oriented normal one-form to H is given by 

r a = L + Cn a (55) 

and so straightforward calculations demonstrate that the associated extrinsic curvature (again in coordinate form) is 

K^dz'dz' = {e7e^W 7 T S )dz'dz' (56) 
= (2Ck v -C v C) dv 2 + (k^l + Cfc^) d6 A d6 B + (2CCj a - d A C) {dvd6 A + d9 A dv) , 

where 



k v = -V a n b V a t 



(57) 
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is again the gauge-dependent quantity that measures how the scaling of the null vectors changes as one moves between 
S v . Since we wish to allow for all values of C, we have not unit-normalized r a . If H is spacelike, the usual extrinsic 
curvature of Section |IT1 defined relative to the timclikc unit normal f a , is 



K 



1 



-K 



(r) 



V2C ' tJ 

Thus, it is clear that the intrinsic and extrinsic geometry of H is fully specified by the set of fields: 



(C, V ,q AB ,K V ,u} A ,k 



AB'^AB) i 



(58) 



(59) 



and so, at least if H is spacelike, we would expect to be able to use these as initial data for evolution into a full 
spacetime. That said, just as qij and Kij are related to each other and the ambient geometry by constraints, our new 
fields are also not all independent. Most obviously, from (38 1 we know that 



£vQab = k 



AB 



Ck 



(n) 
AB 



(60) 



Thus can always be found from C, k£ B and q AB and doesn't need to be independently specified. 

The other relations should be equivalent to the Hamiltonian and momentum constraints: the easiest way to identify 
how is to match equations by the components of the Einstein tensor that appear in them. Doing this we find that 
(|45l-Cx(|46|) gives us 



G ab T a r 



k v ( v) - £ vfyr) + e {n) c v c - d 2 c + 2d A (ccj A ) 

i 



(61) 



~ C \W{n)\\ + 



(n-1) 



9 2 



>) 



which for spacelike H is equivalent to the Hamiltonian constraint. The momentum constraint comes in two pieces. 
First (|45T)+C'x(|46|) gives us 



a_6 



G ab V a r 



£v#(v) 



C 2 \\a 



(n)| 



£ v C-d 2 C + 2ui A d A C 



(62) 



while from (47) we obtain 



(n-1) 



e A G a bT b = C v lj a + 6(y)U) A - d A K V + d B k ( ^l - d A 6 w - Cd A 9( n ) 



(63) 



Subject to these constraints, in future sections we will view (C, q A B, ^ab> ^ab) as an (overdetermined) 

initial data set on H and then use them to perturbatively construct the nearby spacetime. 



2. Deforming H 

Our next step is to understand how the geometry of H changes if it is deformed. As just seen, the intrinsic 
geometry of H is specified by (C,V A ,q AB ) while the extrinsic geometry is also specified if we have knowledge of 
(«y, oj A , k AB , k A g). Already from our earlier calculations we know how to find the deformations of (Jab, k AB , 
and Cj a . Thus we just need to calculate the deformations of C, V A and Ky. Once again we restrict out attention to 
deformations that are normal to the S v and so are of the form: 

X a = at - 7 n a . (64) 

For definiteness we will also need to fix the scaling gauge for the null vectors. We do this by tying that scaling 
to the foliation of H and, once again, it is easiest to see how this works by considering a parameterization of H: 
x a = y a (v,9 A ). Then, as shown in FIG. [2] X a infinitesimally deforms the original surface a coordinate distance Ap 
via 



y a ( v , e A ) -> y a ( v , e A ) + (A P )x a (v, e A ) , 



(65) 
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transferring the coordinates along with the surface. In coordinate terms 

x= w (66) 

The surface coordinates are Lie-dragged by X so we have 

[X, e A ] = and [X, t]=0, (67) 

where, as before, 

ya = _ + g a _ (68) 

The shift vector V A is purely gauge (depending only on the choice of coordinate system on the surface) and so we 
usually simplify our calculations by choosing it to vanish on H itself. However even if we do this, deformations will 
generally cause it to become non-zero off of H. From 

[X,r]=0 => e a A S x V A = -S x V a (69) 



it follows that £ a SxV a = n a 8xV a — and so one can apply (40) to demonstrate that 

6 X V B = d B 1 + (Cd B a - ad B C) - 2(7 - aC)Cb B . (70) 

This rate of change will usually be non-zero. 

The same set of calculations also give us Kx and the deformation of C: 

kx = C v« + a«v (71) 

and 

S X C = C v1 + Ck x -ik v . (72) 

We still need to find the variation of Ky. This is most easily found by expanding the "commutator" Sx^v ~ 5ynx 
using the regular tools. The result is: 

5xKy = Synx + (d B C)(dB(x) + 2£j B (—ds^f + adsC — Cdsct) 
+( 7 - aC) (SQ B u B + n abcd tn b n c l d ) . 

In terms of surface quantities and the Einstein tensor, the Riemann term can be rewritten as 

TP Qa-^b^ccd _ R | ( n ~^) fl a _AB r< Oa^b ^ r < ~ab 

l<-abcd<- U U I — — + - j^ PfflP(n) _ a AB a (n) ~ ^ab^ n - -t^abQ ■ 

We now have formulae for calculating the deformation of all quantities defining the intrinsic and extrinsic geometry 
of if. 



C. Spacetime near H 



With these results in hand we can perturbatively construct spacetime close to H. We will do this using inward- 
oriented Gaussian null coordinates and use the deformation results from the previous subsections to find the perturbed 
metric components. As a concrete example, this construction is implemented for the Kerr-newman spacetime in 
Appendix [A] 

We proceed as follows. Start with an H that is foliated into S v by a coordinate system {v, 8 A } which has been 
chosen so that V A — 0. Scale the null normals to the S v so that V = £ — Cn. This fixes all of our gauge freedoms 
on H . The coordinate system is then extended off of H using the inward null geodesies that cross H tangent to n a . 
Assume an affine parameterization p with the initial scaling set by 



(73) 
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This p will be our oS-H coordinate. Note that these coordinates, though clearly similar in spirit, are not identical with 
the standard Eddington-Finkelstein coordinates familiar from the Kerr-Newman family of spacetimes. As is shown in 
Appendix [A] the geodesies in the standard system are not normal to the surfaces of constant v on the horizon. 

Next, setting p = on H we Lie-drag the coordinates {v,8 A } along the null geodesies to its other level surfaces. 
Thus our full (n+l)-dimensional set of coordinates is {v, p, 9 A }. Because p is an affine geodesic parameter it follows 
that our initial relations: 



n- V 



- 1 and n ■ ea = 



(74) 



are conserved, though in general for p ^ we lose the initial orthogonality between 'f and the e/\. The spacetime 
metric in this coordinate system is then determined by a scalar function C p , an (n— l)-dimcnsional vector function 
Vp and an (n— 1) x (n— l)-dimensional metric tensor function q^ B : 

ds 2 = g al3 dx a dx p = -2dvdp + (2C P + q p AB VpVp)dv 2 + 2q p AB Vfdvd6 B + q p AB dd A d6 B . 

The subscript /superscript ps are included to differentiate these functions from those defined only as initial data on 
H: 



Co = C , q° AB = q AB and V A = V A = . 
These metric-determining functions can be expanded as Taylor series in p. Thus 



(75) 



C 



C + pC' + F -C" 



3! 



<7ab = <7ab+P<7ab 



2| Iab 



C" . . . 

P -5"' 
3! qAB ■ 



(76) 



P V 



/A 



n 2 - 
P_ V „A 

2! 



3! 



V" 



where primes indicate deformations in the n direction (evaluated on H). For example, 

since = there is no zeroth order term in V p . Then, to second order in p, the metric takes the form 

ds 2 = {-2dvdp + 2Cdv 2 + q AB d9 A d6 B } 

+p ^2C'dv 2 + 2q AB V /B dvd8 A + q' AB d9 A d0 B } 



(77) 



(78) 



2 



{2 (C" + g AB V' A V' B ) dv 2 + 2 (g AB V" B + q' AB V B ) dvd6 A + q^ B dd A dd B } . 



We apply our earlier results to calculate each term of the expansions. 

With I = 11 we have a = and 7 = — 1 and so we are dealing with a particularly simple deformation. First, by 
(71) this means that n n = 0. Then (72), (39) and (70) respectively imply that: 



c 



2k A n B and 



V A = 2w\ 



Next ( 73 ) , ( 44 ) and ( 47 ) can be used to obtain the second order derivatives 

C" -- 



3Co A Co A + 1 + e [e) e (n) ~ k A %k AB - \n al3 q ap 



^AB — ""AC B t -A n e B Ly Qy97<5 



1 



-q AB R 1 s'n?n s and 



V 



11 A 



1 / AB iA 

= d B k (n) - d 



(n-1) 
2fc AB ^ B - e Aa Tl a ^ . 



(79) 
(80) 
(81) 



(82) 
(83) 
(84) 
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Often it will be most convenient to leave the expanded metric in the form ( 78 ) . However we can also combine our 
results to write it explicitly as 



ds 2 « {~2dvdp + 2Cdv 2 + q AB d9 A d9 B } + 2p {n v dv 2 + 2Q A dvd9 A + d6 A d0 B } 



(85) 



r (§* A 



OJ U>B 

( n ) B A a. . I 07„(™),~,B 



(n) - k ( h %k™ - \K afj q afi - n a0 £ a n b ) dv 2 ) 



2 (d B 4 - d A 6 (n) + 2k ( ^Cj B - ui A 9 (n) - e%n a0 nP) dvd9 A 
. + ( k Ac k B )C - e%n^eln s C aPlS - j^q^^n 5 ) d9 A d9 B 



This is the metric that we will use in future sections where we construct spacetimes close to horizons. 



IV. HORIZONS 



Before constructing near-horizon spacetimes, we recall some horizon definitions and properties. There are two 
principal types: geometric and causal. Geometric horizons are quasilocally defined and include trapping horizons, 
isolated horizons and dynamical horizons. Like apparent horizons these are defined in terms of the geometry of (n— 1)- 
surfaces and are closely related to trapped surfaces (for which both #(„) < and 9^ < 0). In contrast, event horizons 
are defined relative to the causal structure of the full spacetime. This makes them in one way simpler than geometric 
horizons: one can identify them with just a knowledge of how to calculate null geodesies. However in another way they 
are significantly more complicated: they are highly non-local and defined by the future behaviour of null geodesies. 
For more details and a discussion of the properties of these objects see review articles such as [13 [55J [6"3"HB"6"] or one 
of the original sources as cited below. All of these results are also discussed in some detail in [22] in the same style 
that we use in this paper. 



Geometric horizons 



General cases 



In an (n+l)-dimensional spacetime (M, g a b, V a ), a, future outer trapping horizon (or FOTH) is a n-surface H that 
is foliated by spacelike (n— l)-surfaces (S v , q~AB, eU) such that on each surface: i) 9^ = 0, ii) #(„) < and iii) there is 
a positive function (3 such that 8p n 9^ < [|. These conditions are intended to (locally) mimic those used to define 
apparent horizons [52]: each slice of a FOTH is marginally outer trapped {6if\ = 0) and the other two conditions 
guarantee that it is possible to deform the S v inwards so that they become fully trapped. 

As in our construction of section |IIIB| we can define an evolution vector field 



V a = t - Cr 



(86) 



for each S v , it 



that is both normal to the leaves of the foliation and maps them into each other. Then with 9^ 
follows that £y#(£) = 0. Using (45) this may be expanded into a second order differential equation for C over each 

T,9( t ) < (for some (3) may 



S v . The assumptions that the dominant energy condition holds on the horizon and that 
be used in a maximum principle argument to show that C > 0. In turn this means 

£ v e = -C0 (n) £>O, 



(87) 



since 9< n ) < 0. Under these conditions the FOTH is spacelike or null and the area is non-decreasing. This is the 
second law of FOTH mechanics [2] . 

If C = then the FOTH is null and a type of isolated horizon [35J [35] • Some of the properties of isolated horizons 



will be discussed in section V A Here we simply note that the intrinsic and extrinsic geometries of isolated horizons 
are unchanging in time and there is no flux of stress-energy or gravitational waves across their horizons: they are 
equilibrium states. By the zeroth law of isolated horizon mechanics Kg is constant over an isolated horizon. All Killing 
horizons are examples of isolated horizons. 

If C > then a FOTH is a dynamical horizonjl] and, as noted, is spacelike and expands in area. As for isolated 
horizons, we only summarize properties that are relevant to the current discussion. By another maximum principle 
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argument, it can be demonstrated that the foliations of a dynamical horizon are unique: there is only one foliation 
for which 9^ = [32] . This is very convenient for our discussions as we do not need to worry about whether or not 
geometric properties are foliation dependent. They are, but since the foliation is unique this is fine! On a closely 
related note it can be shown that if a FOTH transitions from being isolated to dynamical it does so all at once. That 
is, on each S v , C is either zero everywhere or it is zero nowhere |32j . 

Dynamical FOTHs (like dynamical apparent horizons) are not uniquely defined. Though the S v cannot be deformed 
within H, they can be deformed out of H (see [22J [32] for theoretical discussions or [3D] for concrete demonstration 
within the Vaidya spacetime). At a local level this non-uniqueness manifests itself in the fact that, for a given S Vl 



changing the scaling of the null vectors will result in ( 45 1 solving for a different V a which in turn will evolve that 



(n— l)-surface into an H' ^ H. By contrast isolated FOTHs are rigid: for that case it can be shown that the only 
allowed deformations are those that change the foliation of H but do not otherwise affect it [H| [33] . 

For the rest of this paper we will consider FOTHs that satisfy the dominant energy condition and so are null 
or spacelike. We respectively refer to them as isolated or dynamical FOTHs. There are also more exotic forms of 
non-FOTH trapping horizons associated with apparent horizon "jumps" (see for example [31] [3j2 [39] ) , inner horizons 
or white holes [5] . However, they can be left aside for the purposes of this article. 

2. Slowly evolving horizons 

With isolated horizons characterizing equilibrium states we can turn our attention to the near-equilibrium regime. 
Intuitively these should be nearly isolated and so nearly null with quantities on the horizon changing slowly in time. 
However, given that dynamical horizon are naturally spacelike, there is no real notion of time on the horizon and 
so the trick is invariantly characterizing this intuition. We do this in the definition below by defining a "slowness" 
parameter e and then using it as a basis of comparison for tracking rates of evolution up H . The characterization 
of a slowly evolving horizon given below is simplified but also slightly strengthened from that originally developed in 
[221 [31] and recently reviewed in [55]. More motivation for the definition can be found in those references. 

Definition: Let AH = {U V S V : V\ < v < V2} be a section of a FOTH with evolution vector field V a — l a — Cn a . 
Define an evolution parameter e via 



e 2 /R 2 H = Maximum C (||c7 (n) f + TZ ab n a n b + df n) /2 



where Rh is the characteristic length scale for the problem. If e <C 1 and the foliation parameter has been chosen so 
that || V|| = \[2G < e, then we say that AH is a slowly evolving horizon (SEH) if on each S v : 

(a) the dominant energy condition holds, 

(b) \R\,Q A u> A , \d A ti A \ and K ab q ab < 1/R%, 

(c) derivatives of a horizon field tangent to the S v are at most of the same magnitude as the maximum of the original 
field. For example, \\d A 9 {n) \\ < 9^ X /R H and 

(d) derivatives of a horizon field "up" the horizon in the V a direction are an order of magnitude e 2 /Rn smaller than 
the maximum of the original field. For example \C v k v \ < (e 2 /R H )n^ ax and |£ V C*| < (e 2 /R H )C max . 

The notation X < Y means that X < k Y for some constant k a of order one while the superscript max indicates the 
largest absolute value quantity over S v for the quantity to which it is attached. 

Let us consider the definition in a bit more detail. First, for standard black holes the characteristic scale is the areal 
radius of the horizon while for black brane spacetimes in an AdS background, it is the radius of curvature of that 
background. Next, note that the definition of the evolution parameter is scaling independent. Geometrically it has 
the implication that the rate of change of the area (or volume depending on the dimension) element is small relative 
to changes in proper length measured up the horizon: 

C 9 e<-^9 (n) e. (89) 
If this invariant condition is met, then by equation (p53| the foliation labelling may always be chosen so that ||V|| < e 



and so that requirement is not independent of condition (88 1. Among other things, this choice ensures that in any 



approach to isolation, C approaches zero in the expected way. Turning to the remaining clauses the energy condition 
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(Ja| has the usual physical implications for the matter fields while (JbJ) and |c]) restrict the geometry on the surface to 
be not too extreme: extreme conditions will generally mean that the horizon will not remain slowly evolving for long. 

Finally Q demands that geometric properties of the horizon change slowly relative to V: this is the clause modified 
from earlier definitions. For simplicity we have opted for a general statement of principal here rather than a list of 
specific quantities that must meet this definition. We have also strengthened the statement by requiring that "time- 
derivatives" of quantities be order e 2 smaller than the original quantities rather than order e. This stronger statement 
continues to be consistent with known examples of slowly evolving horizons. It can also be reasonab ly a rgued that e 2 



(rather than e) is the true scale of slowness for the problem rather than . It is e 2 that appears in ( 88 ) and the rate 
of change of the area/volume element is of order e 2 rather than e. This strengthened statement is required for our 
upcoming demonstration that there is an event horizon candidate close to any SEH. 

An important physical property of SEHs, which supports their interpretation as near-equilibrium states, is that 
they obey versions of the zeroth and first laws of black hole mechanics 22, 34 . For a near-equilibrium state, one would 



expect the surface gravity Ky to be approximately constant and indeed this follows from equation (47). Variations 



WdAKvW ~ (90) 

(thanks to the strengthening of ([d} this is also slightly stronger than the equivalent result in [22l [34] ) . One can also 
combine (|45|) and (46) to derive a first law: 



e<!V^+ r ^r>> (91) 



8ttG J Sv { 8nG 

where we have used the Einstein equation to turn Ricci terms into stress-energy ones [76 . 

B. Causal horizons 

1. General case 

The alternative to geometric horizons are causally defined event horizons. It is well known that these are teleological 
- their position depends on future events. This follows directly from their definition and is probably most easily 
understood with the help of FIG. [3] An event horizon is the boundary of a causal black hole: a region of spacetime 
from which no causal signal can escape. Such a surface is necessarily null and, for outside observers, it is the boundary 
between the unobservable events inside the black hole and those outside that can be seen. One determines the extent 
(or existence) of a black hole by tracing all causal paths "until the end of time" and then retroactively identifying any 
black hole region. The exterior of the black hole is the set of all points for which at least one causal signal reaches J 1 + 
(future null infinity) while the interior is the set of all points from which no signal escapes. The boundary between 
the two regions is the event horizon and it is necessarily a congruence of null geodesies. 

In practice, of course, one cannot trace the paths of all possible null geodesies (let alone all causal curves). However, 
in the case where a spacetime ultimately settles down to an equilibrium state, there is a short-cut to finding the event 
horizon |67j . By the uniqueness theorems the only (3 + 1) dimensional, asymptotically flat, stationary and vacuum 
black hole spacetimes are members of the Kerr family. For these black holes, the location of the event horizon is 
well-known and so once that equilibrium state is reached, one can trace its evolution back into the past to find its 
location at all times (FIG. [3] again). 

Given the nature of their definition, it is perhaps no surprise that event horizons have some unusual properties. A 
particularly important one is that infalling matter can curtail rather than drive the expansion of an event horizon (also 
demonstrated in FIG. [3]). If one thinks of an event horizon as a standard object this behaviour seems counter-intuitive, 
however it is straightforward to see how this happens. 



The event horizon is a null surface and so is ruled by null geodesies. In this case ( 45 ) reduces to the Raychaudhuri 
equation and can be rewritten as 

K t B m - C e 9 (e) = X -d\ f) + \\a m || 2 + K ab £ a i b ■ (92) 

Assuming the Einstein equations and null energy condition the right-hand side of this equation is non-negative. The 
characteristic evolution of the horizon is then determined by the relative magnitudes of KiOm and CiO^y If Cg,6tg\ 
dominates, then 6m necessarily decreases with time and, in particular, infalling matter further curtails the rate of 
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FIG. 3: A schematic that plots both the (spherically symmetric) FOTH and event horizon for a typical Vaidya spacetime in 
which a shell of dust (the shaded gray region) falls into a pre-existing black hole. In this figure, horizontal position records the 
areal radius of the associated spherical shell while the direction of increasing time is roughly vertical outside the event horizon but 
tipping horizontal-and-to-the-left inside. On both sides, inward-moving null geodesies are horizontal while "outward-pointing" 
null geodesies are represented by gray dashed lines. 



expansion. Though at first thought this might seem strange, on second thought it makes sense: the equation is 
simply telling us that the gravitational influence of more mass inside the horizon decreases the rate of expansion of 
the horizon, just as it would for any other set of outward moving geodesies. 

However, in a regime where Kg Om dominates then we have a more naively intuitive situation. In that case an 
increase in the right-hand side of (45) will drive a corresponding increase in K^Om. In particular, if Kg is constant 
or non-increasing, then it is the rate of expansion 6^ that must increase. This behaviour is seen in perturbative 
calculations |69j and we now demonstrate that it is to be expected for near-equilibrium event horizons. 



2. Slowly evolving null surfaces 



Definition: Let AH be a section of an n-dimensional null surface with tangent vector field £ a with characteristic 
scale Rh- Then we say that AH is a slowly evolving null surface (SENS) if for some small e 2 : 

(a) \e\ t) <e 2 (\\a e \\ 2 +TZ ab £ a £ b ) and 

(b) i a can be scaled so that m is of order 1 /Rh and 

Ct9 W~(j^) M- (93) 



(In order to ensure future compatibility with our discussion of SEHs we have written our expansion parameter as e 2 
rather than e). 

The application of these conditions is straightforward. First it is straightforward to see that scaling invariant 
condition (|a|) means that the 6%, term in the right-hand side of (92 1 can be neglected while by (fb]) we can scale the 



null vectors so that on the left-hand the C^9m term may be also dropped. Then we have 



Ke9 w ^\\a {() \\ 2 +TZ ab e a £ b (94) 



and if K£ > 0, the expansion is positive and driven by the flux terms on the right-hand side. 
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These are exactly the kind of conditions that one might expect to hold in a perturbative near-equilibrium regime. If 
the rate of expansion is small then its square will be even smaller while in perturbative calculations one often assumes 
a derivative expansion where time derivatives of quantities are always much smaller than the quantities themselves. 
This is the same kind of assumption as we made in our definition of slowly evolving horizons. 



V. EXTREMAL ISOLATED HORIZONS AND NEAR HORIZON SPACETIMES 



With the mathematical formalism set up and definitions established we are now ready to consider applications. 
The first will be a study of spacetime around an extremal isolated horizon. We begin by considering the phase space 
of isolated horizons. 



A. Phase space of isolated horizons 



Following [35] . an isolated horizon is an n-dimensional null surface with null tangent vector l % and induced (null) 
metric Qij. Its intrinsic geometry is time invariant and so 



£ iQab = 



k W -0 

'MB — u 



9 {e) = and a% = . 



(95) 



However its extrinsic geom etry is als o in variant which give us further constraints on possible on-horizon data. From 



Cik^g = it follows from (42] and (|45j that 



K ab tZ b = and e a A l b e°£ d C abcd = . 



(96) 



Alternatively the restriction on the Ricci tensor follows from (61) and (|62|) - the Hamiltonian constraint plus part of 



the momentum constraint. The rest of the momentum constraint (63) reduces to the zeroth law 



d A K( = 



and a further constraint on the (n+l)-Ricci tensor 



e a A R ab t = . 



Finally in most cases, k AB is fully determined by other quantities. From C ek^g = and (|44|) we find that 



which for Kg ^ can be solved as: 



1 ~ 
2 



RAB = T.^B^ab + d( A Ld B ) + U A U} B 



K AB — 



1 /l 



? a A e b B K ab + d( A uj B ) + luaub 



(97) 



(98) 



(99) 



(100) 



These are our constraints on the phase space of possible isolated horizons - that is the allowed values of 

{qAB, k^l^A, K£, Rab, C abcd } On H. 

The fit between isolated horizon constraints and the Hamiltonian and momentum constraints for spacelike surfaces 
is not perfect, but this is not surprising as we are now dealing with a somewhat specialized null surface rather than 
a general spacelike surface. In particular, data on a single null surface will never be sufficient to fully specify even a 
neighbourhood of that surface: the domain of dependence of a null surface is always empty. 

To get a proper initial value formulation one must specify data on two n-dimensional null surfaces that intersect 
along an (n— l)-dimensional surface [371 EH] ■ For our example one would need: k AB , k^ AB , uj a on an S v along with 
e a A t b e c c £ a l C a bcd orL H and e A n b e c c n d C aoc d on an inward-moving null surface that intersects H along S v (for example a 
congruence of inward moving null geodesies). 

So, in general the data specified on the isolated horizon do not fully specify the behaviour of spacetime, even in a 
restricted neighbourhood. That said if we assume that our expansion around the horizon converges and work at small 
enough p that higher order terms can be neglected, we can approximate that spacetime. We now do this. 
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B. Extremal horizons and near- horizon spacetimes 



There are several equivalent ways to characterize extremal isolated horizons 19J but for our purposes it is most 
convenient to define them as the subset of isolated horizons for which the surface gravity vanishes: ki = 0. Then (97) 
is trivially satisfied while ( 99 1 reduces to 



\rab 



e a A e b B TZ a b + d {A Co B) + u A u B , 



(101) 



(n) 

so that the extrinsic curvature k AB decouples from the constraints and becomes freely specifiable data. 

Applying 9m — K£ — and a AB = to ( 85 1 the spacetime near a extremal isolated horizon takes the form: 

ds 2 w | - 2dvdp + q AB d9 A d9 B } 

+2p ^2Q A dvd9 A + k A l£d9 A dd B } 

( (2w A ^ A + d A u A - K af3 e a n b ) dv 2 

d B k^ B - d A 6( n ) - &A0(n) - e^TZ al3 n^ dvd6 A 
. + ( k AC k B )C ~ e%n^eln s C aPl5 - j^q AB n iS n'n^ d9 A d9 B 



(102) 



-2 (2Jfe^ ) w B 



where the trace of (101) has been used to simplify the dv 2 term that is proportional to p 2 . If we retain only the 
leading order terms in each coefficient this becomes: 

ds 2 w p 2 (2w A w A + d A Cj A - TZ afj £ a n b ) dv 2 - Idvdp + 4pQ A dvd9 A + q AB d9 A d9 B . (103) 

For many readers, this may be a familiar expression. There is a large literature (see |56H58j and references therein) 
on extremal black holes and their near-horizon properties; this metric appears frequently in that work. Very briefly 
it arises in the following way. Start with an extremal black hole in a stationary (n + l)-dimensional spacetime. In a 
neighbourhood of the horizon it is posited [56] that the spacetime metric can always be written in the form: 



ds 2 = r 2 F{r, 9 B )dv 2 + Idvdp + 2rh B {r, 9 A )drd9 B + q BC (r, 9 A )d9 B d9 c 



(104) 



Because the metric is stationary there is no v dependence in any of the terms and the "radial" coordinate r measures 
(affine) coordinate distance from the horizon along null geodesies. The r 2 in front of the dv 2 identifies the horizon as 
extremal. Making the coordinate transformation 



er and v 



(105) 



the near-horizon limit consists of sending e — > 0. Then the p-dependence drops out of the metric and (104) becomes 
ds 2 = p 2 F{0, 9 A )dv 2 + 2dvdr + 2ph B (0, 9 A )dvd9 B + q BC {0, 9 A )d9 B d9 c . (106) 
Viewing this limit as a spacetime in its own right, the Einstein equations reduce to 



R, 



and 



AB 



F 



-h A }iR — dt A h\ 



-h A h A 



(A"B) 



Aq, 



AB 



-y A h A + a. 



Then the equivalence of ( 104 ) with our leading-order spacetime ( 103 1 is obvious. For 

r = — p and h A — —2uj A 



(107) 



(108) 



(109) 



the forms of the metrics match exactly, as do their defining constraint equations ( 101[ ) and (107). 

To leading order the spacetime near a general extremal isolated horizon (including one that might be embedded 
in a non-stationary spacetime) is the same as that near an extremal black hole which is embedded in a stationary 
spacetime. 
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VI. EVENT HORIZON CANDIDATES NEAR SLOWLY EVOLVING HORIZONS 



Next we show that there is always a slowly evolving null surface in close proximity to any slowly evolving horizon. 
In the case where the geometric horizon remains slowly evolving for the rest of time, this null surface is the event 
horizon. To demonstrate this we will use our earlier results to construct the spacetime close to SEH and then search 
for an event horizon candidate in that neighbourhood. 



A. Locating the event horizon candidate 

As for isolated horizons, we base our construction on inward-moving null geodesies with the (v, p, 9 A ) coordinate 
system and work with the metric ( 78 1 . C sets our scale-of-smallness and we consider surfaces £ that can be defined 
by a series of the form: 

P ^ ^p {1) (v,e A ) + P{2) (v,e A ) + P{3) (v,e A ) + P(4) («, e A ) . . . (no) 

where pu) ~ C J and we assume that (similar to other SEH quantities and subject to the inclusion of appropriate 
powers of scaling factor Rh ) 

^<C J+1 and \\d A p {J) \\<C J . (Ill) 

Then, to second order, the induced metric on £ is 

dS 2 = q^ B d6 A d9 B (112) 
+ {2 (C + P{1) C') dv 2 + 2 (p (1) q AB V B - d A P(i)) dvd0 A + (p ( i)&„) d6 A d0 B } 

+ 1 Up {2) C' - 2 P{1) + p 2 {1) C" + p\ x) qwV lA vAdv 2 



IB 1 2 -/ i~i/B 1 1 „2 ~ <V/B 



+2 [ -d A p {2) + p (2) <?abV + -p (1) q AB V + ~p (1) <ZabV 
P(2)9ab + ^i)«Ab) d0 A d6 B 1 



where = £vP(i) = dp^/dv. 

Any event horizon candidate will be null and so we need to solve for the pi n so that the determinant of this induced 
metric vanishes. To that end recall that a general metric of the form 

dY? = Fdv 2 + 2V A dvd6 A + h AB d8 A d8 B (113) 

has determinant 

(F - /i AB W 6 ) x det(/i) . (114) 
If h AB is spacelike (as <jab is in our case) then the full metric determinant vanishes if and only if 

F-h AB V A V B = 0. (115) 



We can apply this to the induced metric (112) and solve order-by-order for the pij)- To zeroth order any such £ is 
already null but at first order we must have 

C + P(l) C' = {) (116) 

and so find that 

P(i)=-— ■ (H7) 
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The determinant also vanishes to second order if 



2p (2) C" - 2p (1) + p^C" + p^abV^V 
This is easily solved for p( 2 ) 

P( 2 ) - ^7 I P(i) - §P(i)C" - P(i)V A d A P(i) - - || dp (1) 



(118) 



(119) 



To get this in terms of geometric quantities we substitute in ( 117) and apply the zeroth law for SEHs. Thus to second 
order 



C 



+ \ \Ck v - Ck v ~ \c 2 C" - 2Cu A d A C+ \\\dC\\ 2 ) , 



Ky «y 



where from (82): 



C" 



-3(D A (I'a ■ 



+ y _ AB (n) _ 2 ~ Q/? ' 



(120) 



(121) 



This is a null surface and in fact this is the only null surface that lives entirely in the regime of the near-horizon 
approximation. There are other null surfaces that pass through the region but this is the only one that remains there. 
It is clear that in the isolated limit — > and so one can reasonably argue that if the horizon remains slowly 
evolving for its entire future and it ultimately asymptotes to isolation, then this is the event horizon. 



B. Properties of the event horizon candidate 



With our accumulated computational infrastructure it is straightforward to find the geometrical properties of this 
near-SEH null surface and demonstrate that it is a slowly evolving null surface. We could work directly from the 
metric and surface defined by ( 120 1 but it will be more convenient to calculate the required terms as deformations 
generated by the vector field: 



X 



c_ 



First the induced metric and area/volume element on the S v that foliate £ are 



-(f) 
9ab 



r(£) 



<Zab 



— I &ar and 



c_ 



*(n) 



(122) 

(123) 
(124) 



As would be expected, to leading order these match the corresponding quantities on H. This is also the case for other 
geometric properties with the exception of the expansion which is already at sub-leading order. In that case 



<?(V) + Sx&(V) 
-C0(„) + C0( n ) 



(125) 



(h (e) \\ 2 + G ab £ a £ b ) 



(Ik 



\ 2 + G ab e a £ b ) 



The first line is just the leading order expansion while the transition to the second applies (45), (46) and (72) and 



properties of slowly evolving horizons to identify and discard higher order terms. The transition to the final line is 
then obvious. In any case we recover a "first" law for £ : 



A£)n{£) 



+ G ab tt 



(126) 
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I I (£\ 

where we have also applied (73 1 to demonstrate that to leading order k\ m Ky- This is a direct demonstration of 
how the first law holds on £. However it would also be straightforward to apply the deformations to check that this 
is a slowly evolving null surface and so must obey a first law. 

Note however, that even at leading order the expansion of £ differs from that of H . Combining (45 1 + C x (46 ) and 
applying the slowly evolving conditions it follows that: 

n V 0(v) ^d B (d B C-2CQ B ) + \\a (e) \\ 2 + G ab £ a £ b . (127) 

That is 

M(V) ~ K e £)e (fj ~ d B {d B C - 2Cu B ) . (128) 

At leading order the expansions differ by a total derivative. The origin of this term can be better understood starting 
from our defining relation: [f , X] — 0. Expanding this gives 

[V + V,X] = => 6 X (6 V 1) = S v {S x e) - 8 {x tfi. (129) 

Thus (again applying properties of slowly evolving horizons) 

5 x (e9 v ) ^d B (S X V B ) (130) 



which leads to ( 128 ). The total derivative term is the divergence of the induced shift on £. If C is constant (or nearly 
constant) over a slowly rotating surface, the expansions match at leading order. However in general, from the point 
of view of expansions, the rotation and/or non-constant C induces a mismatch in points on H and S v . In some sense 
this is an issue of choice of coordinates: if we integrate the expressions over the S v and so compare rates of change of 
area, then the total derivative integrates out (for closed S v ) and they do match. 

C. Examples 

Event horizon candidates have been seen previously for certain spacetimes. We now compare our general result 
with those specific ones. 

1. Event horizon candidates in Vaidya spacetimes 
We begin with VaidyapTJ [54] spacetimes. The Vaidya metric 

ds 2 = - M - '^y^j dv 2 + 2dvdr + r 2 (d6 2 + sin 2 6d(f) 2 ) (131) 
describes a spherically symmetric black hole that is being irradiated by infalling null dust with stress-energy tensor 

T ab = ^NaNt • (132) 

Any non-decreasing mass function will satisfy the energy conditions. Outward and inward-oriented null vectors are 
given by 



d_ 
dr 



(134) 



Then for general r 



There is a FOTH located at r = 2m. 



V — 2 777, 2 

9 W = — and 9 {n) =--. (135) 



r 
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We now specialize to the case of a black hole that transitions from an initial mass mi to a final mass m<i — 2m\. 
Then Rh = mi and we scale: 

r = Rmi , v = Vnii and m = Mm\ . (136) 
The horizon is slowly evolving if M <C 1 [HI [53] (here and forthwith the overdot indicates a derivative by V). Now 



from ( 133 ) outward oriented spherically symmetric null surfaces must be solutions of 



dR 
dV 2 



If we assume a hierarchy of derivatives so that M ^> M 7$> M ^> M . . . then we can search for perturbative solutions 
of the form 

R = 2M (l + aM + (3M 2 + 7M + . . . ) . (138) 

Implementing this and switching back to unsealed coordinates, it turns out that to second order there is a null surface 
at 

A £) w 2m + 8mm + 32m(2m 2 + mm) . (139) 

This direct calculation matches our more general one. For this scaling of the null vectors it is straightforward to 
find that on the horizon 

k v = — and C = 2m (140) 
4m 

and since the FOTH at r = 2m is spherically symmetric 

R=-^. (141) 
2m/ 



Further, —r is an affine parameter for the inward moving null geodesies. Setting p = — (r — 2m) we apply (120) and 
find that 

pM\ — —8mm and (142) 
P{2) — —64mm + 32m 2 m . 



Converting back into regular r coordinates, this becomes (139). 



2. Boost-invariant black brane spacetimes 

Our next example is the five-dimensional black brane spacetime that is the fluid-gravity dual to Bjorken flow[53j. 
The spacetime metric takes the form 

ds 2 = -r 2 A(f)df 2 + 2drdr + (1 + rf fe h ^ ' r Uz 2 + r 2 e c(f ' r) (dx 2 + dy 2 ) , (143) 

where the defining functions can be expanded as 

j4 <V,_, W „) + *M + *M + ... (144) 

K f , r) _ W „ + ^ + hW + ... (1 45) 

o(f,r) = cM + C ^ + ^ + ..., (146) 

and t > 1 (so the expansion parameter f~ 2 / 3 « 1). At the same time r is taken to be sufficiently small so that 
v = rf 1 / 3 is always of moderate size. 
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Applying the Einstein equations, one can solve order-by-order for the defining function. At lowest order 

7T 4 A 4 

A (v) = 1 — , b a (v) = and c (v) = (147) 

while at higher orders things become considerably more complicated. Here we will not be concerned with the details 
of these calculations: the results that we need may all simply be read out of |53j . In particular it is shown there that: 

A f 1 1 / 1 \ 1 / 29 5 171og2\ 1 

Teh ^^+^ 1 + ^^r 2 + -j+^^ — ^^j+...J (148) 

A f 1 1 / 1 \ If 25 1 251og2\ 1 

tab = ^(r 0+x ^r 1 + x ^ ^ 2 + _j +x ^^3--_ + _- 1 ^j + ...| (149) 

where 

r = it (150) 

*--i-T < 151 > 

A<5 2 1 log 2 . , 

r 2 = — 5— (152) 

2 3 24 18ttA v 7 

1 A<5 3 log(TrA) C 71og 2 2-121og(7rA) 

r 3 - + -3 + 7^2 + 7^2 ( 153 ) 

for non-trivial constants #i, <5 2 and <5 3 . Then a straightforward subtraction gives 

1 fj_\ L_ M 1 | log2 



r EH -r AH = 7 ^ ___»__ + __+] . (154) 



We can compare this with the predicted horizon separation from current calculations. It is also shown in |53] that 

1 1 /log 2 1 

9^2 ~ Af8 /3 I ig^" ~ 54 



O-X-^i^-^... (155, 



and 



If these are substituted into (1201 then we should recover (154 1. The calculation is fairly straightforward. Because 



the metric is vacuum and has the planar symmetry (119) reduces to 



P(2) = ^3- (ck v - Ck v + -C-a AB a J ;^, 



Ck v -C KV + -C 2 a { ^ ) ) . (157) 



The shear term is of lower order than the others and so may also be neglected (||ct^|| ^ \C for a slowly evolving 
horizon). The only complication arises from the fact that both Ky and C are themselves expressed as series. Thus 
itself contains lower order terms, some of which are of order p( 2 ). Once these are properly accounted for, we recover 
( 154| ) as expected. 



VII. DISCUSSION 



Given geometric data on a (foliated) n-dimensional hypersurface H, the main work of this paper was to pertur- 
batively reconstruct the (n+l)-dimensional spacetime near that surface. We did this in a Gaussian null coordinate 
system determined by the foliation of H . To second order in the affine parameter on the inward-oriented null geodesies, 



the resulting metric was (85). At leading order this is entirely determined by the intrinsic geometry of H while at 
next order the extrinsic geometry also contributes. However, at second and higher orders more knowledge (in the 
form of the Ricci and Weyl tensors and their derivatives on H) is needed. These initial calculations were entirely 
geometric and neither assumed that H was any kind of horizon nor made use any field equations. They apply to any 
hypersurface H (that can be foliated by spacelike surfaces) in any spacetime. 
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Of course if H is spacelike and the vacuum Einstein equations hold, then it is well-known that there is a good 
initial value formulation of general relativity based on the intrinsic and extrinsic metric of H (Section pi). In that case 
the intrinsic metric and extrinsic curvature are sufficient to determine a series expansion of the metric to all orders 
(in timelike geodesic normal coordinates). This appears to contrast with our calculation where, at second order, 
knowledge of the Weyl tensor is also required. However, a little investigaiton shows that the contradiction is only 
apparent. If H is spacclikc (C =/= 0) and we assume that the Einstein equations hold, then the required components 
of the Weyl tensor can be found from the derivative of the inward-oriented extrinsic curvature "up" H via equation 



(44): 



e a A n»e c B n d C abcd = ^ (m£2 + M&b) ( 158 ) 



, (n) ,{n)C _ Q7T - rp c c 
AC B 7 TTQAB^cdii n 

[n— 1) 

1 /l s 1 



C\2 RaB ^2 



h») k AB + 6 {t) k AB 



» Al)C _ t o a c b 
V C(A B) 



~ A e b B 1Z a i ) + d( A Co B ) — Q A u> B ) 



Essentially this is a constraint equation on H that determines e A n c e B n d C ac bd as a function of the intrinsic and 
extrinsic geometry terms. Though we have not done the calculation explicitly, constraints should similarly determine 
the higher order quantities. 



If H is null the situation is different. In that case C — and ( 158 ) is not well-defined. However this is not surprising 
since in this case we would not expect a good initial value formulation: the domain of dependence of a null surface 
is empty. Physically this is because extra information that has travelled "parallel" to H can influence spacctime 
arbitrarily close to H. As mentioned earlier, to get a good initial value formulation, data must be specified on a pair 
of intersecting null surfaces [68 . Thus apart from on an isolated horizon one would also need to specify data on, 
for example, an outgoing past-directed null cone based originating from some S Vo . Specifically in vacuum one would 

need qAB, k AB , and Cj a on S v along with e A £ c e B £ d C ac bd on H (this must vanish for an isolated horizon) and 
e A n c e B n d C a cbd on the past-oriented null-cone [5J[37]. This would be sufficient to specify the spacetime for v > v Q 
and p < (at least while the coordinates are well-defined). 

We considered two applications of these results. The first was a straightforward reconstruction of the spacetime 
near an extremal isolated horizon. At leading order this can be done without any need for off-horizon information 
and the result is a near- horizon spacetime in the sense of [56, 57 . In fact, these near-horizon spacetimes are actually 
exact solutions of the Einstein equations: in this case, throwing away sub-leading corrections turned an approximate 
solution into an exact one! The details of exactly why this happens deserve further consideration but we leave this 
small puzzle for later investigation. 

The second example was more involved and we demonstrated the existence of an event horizon candidate which 
hugs any slowly evolving horizon. This general result was foreshadowed by similar results in the case of Vaidya and 
near-equilibrium black brane spacetimes, however we believe that this is the first time that it has been demonstrated 
in full generality. It should be emphasized that this tentative identification does not violate the teleological nature of 
true event horizons: in order to identify this SEH-hugging null surface as a true event horizon we must assume that 
the trapping horizon remains slowly evolving for the rest of eternity and that it ultimately settles down (or at least 
asymptotes) to equilibrium. If both of these are true, then we have identified a null surface that asymptotes to that 
equilibrium state and so is the event horizon. 

Of course in all of this work we have assumed that the series expansion of the metric converges. We expect this 
to be the case. Higher order derivatives of our geometric quantities will also depend on surface quantities and their 
derivatives and we do not expect them to blow up. However, this is not a completely trivial result. The first attempt 
to demonstrate the existence of an event horizon candidate near slowly evolving horizons used the standard (n+1)- 
formalism from Section [il] to expand spacetime around the (spacelike) dynamical slowly evolving horizon |55j . At 
leading order the results of that work essentially matched those found here. However subsequent calculations have 



demonstrated that the "second" order corrections in the standard expansion are actually proportional to 1/VC: in 



timelike geodesic coordinates the series analogous to (110) does not converge. In retrospect this is not so surprising. 



Relative to our Gaussian normal coordinate system, the inward oriented timelike normal to a slowly evolving horizon 
is 

1 / a \ a Tn / Ft \ a 

(159) 



and so time derivatives relative to this vector pick up factors of (2C) -1 ' 2 . Intuitively small changes in proper time 
along the normal geodesies can correspond to large changes in the horizon coordinate v and therefore significant 
changes in the usual geometric quantities. 
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In conclusion we briefly consider future applications of this formalism. It is easy to see that "just outside" the 
event horizon candidate there will a timelike surface that similarly hugs any slowly evolving horizon. Such a surface 
could be treated as a stretched horizon from membrane paradigm [70] • Via this link much of the membrane paradigm 
formalism will be translatable into the language of slowly evolving horizons and their accompanying surfaces (and 
vice versa). We expect that new astrophysical insights may follow from this cross- fertilization with obvious targets of 
study including black hole ring-downs [72 j. spin-flips |46j post-merger recoils [71] and anti-kicks [38] [43]. It will also be 
possible to write many of the membrane paradigm results in a form that will apply to general black holes and branes 
in general dimensions. Connections with blackfolds [73] are an obvious target of study. 
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Appendix A: Kerr-Newman in EF-normal coordinates 



In this appendix we demonstrate the construction of the Gaussian null coordinate system for a Kerr-Newman 
horizon. To begin, recall that in standard Eddington-Finkelstein coordinates the Kerr-Newman spacetime takes the 
form: 

ds 2 = - (l - *^*\ dv 2 + 2dvdr 2a(x-A)sin 2 ^ # ^ 
-2a tiJOdrd* + ^9- + ^(x 2 -a 2 Asm 2 ^ 2 

where A = r 2 — 2mr + a 2 + Q 2 , \ = r 2 + a 2 and £ = r 2 + a 2 cos 2 0. The isolated horizon H is at r Q , the larger root 
of A(r) and we consider the foliation S v of surfaces of constant v. Note that though based on ingoing null geodesies 
there are two things that distinguish these standard coordinates from Gaussian null coordinates: 

1. On H , the dvdfi-term is non- vanishing. That is, is not orthogonal to the foliation surfaces of constant v and 
as such it is not a V a candidate. 

2. The coordinate r is an affine parameter for a family of ingoing future oriented null geodesies but that family is 
not orthogonal to the S v (as evidenced by the non-zero drd<p term). 

The first of these difficulties is easily resolved by the coordinate transformation 

cP = <p+(~)v (A2) 

with Xo = r o + a 2 which "unwinds" the S v on the horizon. The metric becomes: 

^2 { Z 2 A-a 2 sm 2 e(r-r ) 2 \ ^ 2 2S 2asin 2 g {2 2 , 

ds = - ^ ) dv + ~^ dvdr + S(r 2 + Q 2) ( AE ° + *( r ~ r o)) ( A3 ) 

„ . 2„, ^ ,„2 sin 2 6»(Y 2 - a 2 Asin 2 6>) , ,, 
-2a sin 2 9drdcj) + zZd8 2 H y — — '-d(j> 2 

Zj 

where E = r 2 + a 2 cos 2 #. Then the full (null) three-metric on H is 

dS 2 = Z d9 2 + ^^# 2 (A4) 

which is now of the expected form. 

The second difficulty is resolved by rewriting the metric relative to the correct set of geodesies. Relative to the 
current coordinate system a suitable pair of (cross-normalized) null normals is 

. d , / a 2 sin 2 fA 8 fr 2 +a 2 \ d ( a \ 8 

e= d~v and n = -(^2Y7- )to-(-^r )fr-{*T* )w (A5) 
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Keep in mind that these are defined only on the horizon. Now consider the family of null geodesies that crosses H 
with n as its tangent vector field. We identify them by the point (u, 0, <p) where they cross H and parameterize them 
with affine parameter p so that p — on H and increases inwards. We perturbatively construct the geodesies up to 
third order in p: 



dX c 



dp 



p=0 



p^_ d 2 X a 
2 dp 2 



p=0 



p 3 d 3 X a 
~6~ dp 3 



(A6) 



p=0 



On the right-hand side of this equation and henceforth, the labelling subscript (v, 9, ip) is omitted but understood. 
The first two coefficients are trivial: 



X' 



\H 



[v, r„, 6, ip] 



(A7) 



and 



dX a 
dp 



H 



a 2 sin 2 £ 



2S„ 



,0,— 



(A8) 



The next two follow from the geodesic equation: 



n^V a n a = ^ 



d 2 X a 



dp 2 



H 



(A9) 



and 



n 7 V 7 {n P V B n a ) 







d 3 X a 



dp 3 



H 



(A10) 



The right-hand sides follow from expanding the left-hand side and making appropriate substitutions from the earlier 
derivatives. The Christoffel symbols f 1 " ^ a nd their derivatives dpT a ^ s only need to be evaluated on H . 

Once these quantities are calculated (A6| defines a transformation from (v,r,9,ip) to (v,p,9,tp) coordinates. Then 
the second order expansion of the metric is: 



9ap ~ g a p + pg al3 + y 5 Q/3 



(All) 



where the zeroth order components are 



(0) 



JO) = ( *o 



(A12) 



sm 



the first order corrections are found to be 

A' 



9 (1) 



o (1) 



I oo 

(i) 



Xo 



2a 2 sin9 cos9 



So 

asin 2 6> 



A' 



S 

2a 3 Xo sin 3 0cos0 



(A13) 



2ar Xo sin 
£2 
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and the second order corrections are (the somewhat complicated): 



9 (2) 

■ ' r r 



(2) 
9v9 



9 (2) 



(2) 
998 



(2) 

9 dip 



a 2 sin 2 0\ , 2 (r (2xo + a 2 sin z #) 



4x 2 So 



(A 



,2^2, 



XoS 2 



A' 



S 2 - 4a 2 r 2 sin 2 ( 



(A14) 



a 2 sin# cos#(3x 2 + a 2 sin 2 #) \_ , ( 2a 2 r Q s'm9 cos6>(x 2 + a 2 sm 2( 



a 3 sm 4 6 
4XoS 2 

9 • 2 / 

a r c sm t 
2E? 



2XoS 2 
(A') 2 - 



A' 



r c sin 2 0(4x o + 3a 2 sin 2 6») \ / a Xo sin 2 6»(2E 2 - a 2 sin 2 6> (5r 2 - a 2 cos 2 6»)) ' 

+ (1 + cos 2 6»)a 2 r 4 - (5cos 4 + 7cos 2 6> - l)a 4 r 2 + sin 2 6»cos 2 6>(cos 2 6> - 5) a' 

El 



fl3 (4Xo - S ) sin 3 0cos0"\ , f a 3 r x sin 3 6> cos6»(E + 6a 2 sin 2 ' 



9 (2) = 

tfiplp 



2sg ; v ^4 

a 4 sin 6 6'\ , 2 | f a 2 r x sin 4 6'(3x + 4a 2 sin 2 ' 

^sFj ( } " 



s 4 



A' 



+ 



X 2 sin 2 6» (r° + (5 cos 2 <9 - l)a 2 r 4 + (11 cos 4 (9 - 15 cos 2 6> + 5)a 4 r 2 + sin 2 6 cos 2 6»(l + cos 2 (9)a 6 



These same expressions are also found when calculating directly from (78). This, of course, isn't a surprise but it 
does provide a reassuring cross-check on potential typographical errors. 
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